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Abstract 

We give a matrix formulation of the Hamiltonian structures of constrained 
KP hierarchy. First, we derive from the matrix formulation the Hamiltonian 
structure of the one-constraint KP hierarchy, which was originally obtained by 
Oevel and Strampp. We then generalize the derivation to the multi-constraint 
case and show that the resulting bracket is actually the second Gelfand-Dickey 
bracket associated with the corresponding Lax operator. The matrix formu- 
lation of the Hamiltonian structure of the one-constraint KP hierarchy in the 
form introduced in the study of matrix model is also discussed 
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I. INTRODUCTION 



The Gelfand-Dickey (GD) hierarchy is defined by 



l n = d n + u n „ 1 d n - 1 + 



+ u 



which satisfies the hierarchy equations 

dl n 
dt k 



[{In) 



k/n j 1 
+ j l n\ 



;i.2) 



Here A± denote the differential part and the integral part of the pseudo-differential operator 
A. The second Hamiltonian structure of (|1.1| ) is described by the GD bracket which in 
operator form can be written as 



,5H, 



e^(— ) = {l n ,H} = {u n - 1 ,H}d 

0l n 



,n-l 



+ { Un -2,H}d n - 2 + --- + {u ,H} 



5H 



5H 



[I'll rj S ) + ln lfl\ f-j I 
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where 



$ln 5u 



5H 



6u 



;i.3) 



;i.4) 



n-1 



Note that, from (|1.2p , the equation for u\ is trivial thus we can set u\ = 0. However, 
imposing such a constraint leads to a modification of the GD bracket ( |1.3| ) due to the Dirac 
reduction. The modified bracket is found to be 



- rn .5H. 5H. , , .5H , . l ri r ri 5H, 



1.5) 



where iti in Z n and in 5H/5l n are both set to zero and res(X)i &id l ) = a_i. 

The above description of the GD hierarchy is relied on the use of fractional-power pseudo- 
differential operators P] associated with the scalar Lax operator l n . However, one can put 
the formalism on a more general setting by considering the GD hierarchy as a reduction of 
a system of n first-order equations. Following Dickey ||, one can substitute for the scalar 
Lax operator l n a n x n first order differential operator 



C 



U 



Id + U 
( -1 



V u «1 



d = d/dx 

\ 

-1 



... -1 

U n -l J 



(1.6) 
(1.7) 



This is the basis of the extension to the more general situation when the matrix U belongs to 
a semi-simple Lie algebra, a task that was accomplished by Drinfeld and Sokolov [|J . Using 
matrix notation, the hierarchy equations (|1.2j) can be expressed as 
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dtC=[Q,C] 



(1.8) 



where Q is another n x n matrix which can not be arbitrarily chosen. We have to choose 
properly the matrix Q such that consistent with the form of dtC. In fact, in view 

of ( [Op , if the last column of Q is given then the rest of the elements of Q can be fixed. 
Following this strategy, it has been shown that the Hamiltonian structure of the GD 



hierarchy can be extracted from (|1.8j) and turns out to be the GD bracket ([L 

Recently, there has much interest in the so-called constrained KP hierarchy [§-|13[ which 
has a pseudo-differential Lax operator of the form 

m 

L M =d n + u n ^d n - 1 + ■ • • + u + J2 <Pkd-^ k (1.9) 

k=l 

and satisfies the evolution equation 

i"'"^ — [{L(n,m))+ n , £(n,m)] (1-10) 

Otu 



From ( |1.10|) , it can be shown that each <pi (each ipi) is an eigenfunction (adjoint eigenfunction) 



of the constrained KP hierarchy, i.e. 



f| = «£ M )f>i) H = -m n , m) )T^) (i-ii) 



The bi-Hamiltonian structure of ( |1.10p have been constructed by Oevel and Strampp 



[see also |f§). The matrix formulation of the cKP hierarchy and its Hamiltonian structure 
have been discussed in refs. [OLinj]. In contrast to the affine Lie algebraic approach 
in this paper we shall follow Dickey's approach || to give an elementary derivation of 
the Hamiltonian structure associated with the Lax operator L( ra m ). In our approach, the 
Hamiltonian structures obtained by Oevel and Strampp || come out quite naturally. Since 
the first Hamiltonian structure can be obtained from the second Hamiltonian structure by 
replacing the Lax operator L( n>m ) by L( n>m ) + A, where A is called the spectral parameter, 
we shall focus only on the second structure. 

Our paper is organized as follows: In Sec. II, we follow the approach close to that of 
@ to build up the matrix equation ( |1.8|) associated with the Lax operator Lr n \\ and then 
derive the Hamiltonian structure which was obtained by Oevel and Strampp. In Sec. Ill, we 
generalize this formulation to the Lax operator L( n ,m.) an d obtain its Hamiltonian structure. 
We further show that this Hamiltonian structure is, in fact, the GD bracket ( |1.3| ) defined 
by the Lax operator L( n ,m)- 111 Sec. IV, we work out a few simple examples explicitly. 
Concluding remarks are presented in the Sec. V. 

II. MATRIX FORMULATION OF OEVEL AND STRAMPP'S HAMILTONIAN 

STRUCTURE 

Now we consider the matrix formulation corresponding to the Lax operator for the one- 
constraint KP hierarchy 

= d n + Wn-i^ -1 + ■■• + % + 0<9~V (2-1) 
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u = 



(2.3) 



A matrix representation can be easily found by expressing the constraint equation L( n ^ip = 
in a matrix form as £( n i)$ = 0, where is a square matrix and $ is a column matrix. 

An ansatz is given by 

£ (nl) = Id + U d = d/dx (2.2) 

/ -ip ••• \ 
•• -1 o 

-1 



-1 

V u ui U n -i J 

where U is a (n + 1) x (n + 1) matrix and the numeration of rows and columns will be 
hj = -l,0,l,2,---,ra- 1. 

We introduce another (n + 1) x (n + 1) matrix Q and ask the commutator [Q, £( n ,i)] = ^ 
to be consistent with the form of <9t£( n ,i)- Let's now compute the matrix elements of the 
commutator: 

(1) % = -1, or j = -1 



N-i-i = -Q'_! _! + 0Q-l,n-l + ^Q 0> -1 

= - eiQ-i,*-i + ^Qo,i + «»Q-l,n-l (eo = ^, Cl 

A^,_i = -o; _! + (f)Q l ,n-i + Qi+i,-! (i = 0, 1, 2, • • • , n - 2) 

n— 1 

Nn-1,-1 = -Qn-l -1 + <f>Qn-l,n-l - <PQ-1,-1 ~ U jQj-l 

3=0 

(2) i ^ -1 and j ^ -1 



(2.4) 

e„_ 2 = 1) (2.5) 
(2.6) 

(2.7) 





— — Qi,0 ~~ i^Qi-l + u oQi,n-l + Qi+1,0 


Wn-1,0 


— ~~ Qn-1,0 ~~ i^Qn-l-l + 'WoQn-l.n-l ~ 




Qi,j—1 V'jQijfi—l Qi+l,j 




= ~Qn-l,j ~~ Qn-l,j-l + u jQn-l,n-l ~~ ( 


where i — 0, 1, • • 


• , n — 2 and j = 1, 2, • • • , n — 2. 



n— 1 

-i,o - X] u j'Q?,o 

J=0 



n-1 



fc=0 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



In order to make N of the form d t C( nt i), we write iV = 5£( n ,i)\ i- e - 



iV_ 1>0 = -Sij) 

iVn-1,-1 = 50 



(j = 0, n-1) 



(2.12) 
(2.13) 
(2.14) 



and iVjj = for other In the matrix formulation the equation 6£( n ,i) = [Q, £{n,i)] 

would serve as a definition of the Hamiltonian structure. In other words, 5ip, 5<p and 8uj 
are to be identified as {ip,H}, {<p,H} and {uj, H}, respectively. Of course, certain matrix 
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elements of Q must be identified as the components of the gradient of H. However we shall 
do such an identification later. 

We shall see that one can solve Q it j in terms of Qo-i and Q^ n -\ (i = — 1, 0, • • • , n — 1). 
To this end we introduce 



n-l 



n-l 



Qi = Y, Q^ k 4 = E Sr x ~ h Qhj 

k=0 k=0 

n-l 

N = N ik d k 

k=0 



then (|2.8|) and ( |2.10|) is equivalent to 

Ni = -dQi - ipQi t -\ + Qi+i + Qi, n -i(L(n,i))+ (i = 0, 1, • • • , n - 2) 
While (fl9|) and ( fTTTD give 



n-l 



N n -i = -dQ n -i - i)Q n -\-\ - <PQ-1 + Qn-l,n-l(L( n ,l))+ ~ E U hQk 



(2.15) 
(2.16) 

(2.17) 

(2.18) 



k=0 



Since Ni = (z = 0, 1, • • • , n — 2) we deduce a recursion relation for Qt from ( |2.17|) 

4+i = 0&- Qi,n-i(L {ntl) )+ + ipQi,-i (i = 0,l,---,n-2) (2.19) 

Define 



n-l 



4 = ~ E M ^4 ~ ^n-1 ~ 



fc=0 



Then (|2.19|) also holds for i = n — 1. Hence 

4 = <9Qo — Qo,n-l(^(n,l))+ + V'Qo-l 

4 = #4 - Qi in _i(L (n) i)) + + V'Qi-i 

= <9 2 Q - <9Q ,n-l (£(«,!))+ - Ql,n-l(^(n,l))+ + ^Qo.-l + V'Ql-l 



i-1 



i-1 



4 = 9*(Qo - E ^ 1_fc Qfc,n-l(^(",l))+ + E 9" 1_fe Qfc,-lV') 



fc=0 
n-l 



fc=0 
n-l 



Qn = o n (Q - £ a- 1 - fc g fc , n _ 1 (L (nil) ) + + £ d- 1 -^^) 

k=0 k=0 

By the virtue of ( |2.20| ), Q n must be a differential operator of order < n — 1 it follows 

n— 1 n— 1 

(4 - E ^ 14 ft ! „-i(i ( n,i)) + + E d- 1 - k Q k ^) + 



(2.20) 

(2.21) 
(2.22) 
(2.23) 

(2.24) 



fc=0 



k=0 

n-l 

(Oo — E ^ l fe< 5fe,n-l(-^(n,l)) + )+ = 
fc=0 



(2.25) 



Define 

X = Y, d~ X ~ k Qk,n-i = Q n ~ X (2.26) 

fc=0 

we get 

Q = (X(L M ) + ) + = (XL M )+ (2.27) 

and 

4 = 9 i (XL (n , 1) ) + - (d l X) + (L {nA) ) + + (d'Q- 1 )^ (2.28) 

From ( ]2T20| ) 



n-1 

- £ u k [d\XL M ) + - (d k X) + + (d'Q- 1 )^} 

k=0 

= -(L M )+(XL M ) + + {L {nA) X) + {L (nA) ) + - (L^Q- 1 )^ - 0Q_! (2.29) 
In view of 5uj = N n _ij we have 

<5(L( ni i)) + = (L( ni i)X) + (L( n) i)) + - (L(„ i i)) + (XL(„ ) i)) + - (L( ni i)Q" 1 )+V ; - 4>Q-i (2.30) 



It remains to characterize Q-i and Q Since Ni-\ = (i = — 1, 0, • ■ ■ , n — 2), (|2.4j ) and 
(ED S ive 

= f(<PQ-l,n-l + VQo,-i) (2.31) 
= <5i,-i - <t>Qi,n-i (i = 0,l,-..,n-2) (2.32) 

If we define 

n-l 

<5n,-l = -0(5-1-1 - 51 U k Qk-l ~ iVn-1,-1 (2.33) 

A;=0 



f2.7| ) becomes ( |2.32| ) with z = n — 1. It follows that 



(d- n Q n ,-i + ■■■ + ^Qi.-i) = (Q" 1 )' - Q n ~ V- (2.34) 
Using (Q- 1 )' = [d,Q- 1 ] we get 

Q n ,_i = d n (<5o,-i " <5~ X <9 - ^0) (2.35) 
Since Q n ,-i is a scalar, it follows 

[d n (Qo,-i-Q~ 1 d-X ( J ) )d- 1 } + = (2.36) 
Moreover, since d n Q~ 1 is a pure differential operator we have 



d n Q~ 1 = [d n (Q ^ - Xct>)d- l } + (2.37) 

Next from (]2j) 

Q-i,i-x = -Q'_ M + VQo,» + «iQ-i,n-i (i = l,2,.--,n-l) (2.38) 
Define = —Q'-x o + V'Qo.o + woQ-i,n-i then in view of (|2.38|) we write 

= -(Q-l)' + ^4) + Q-l,„_l(L(„,i)) + (2.39) 

which implies 

Q-i,-i = -dQ-i + ^Qo + Q_i, n -i(L( n ,i)) + (2.40) 

Since a scalar is invariant under the adjoint operation; i.e. / = /*, applying the adjoint 
operation to the both sides of ( |2.40| ) gives 

Q-L-i = (Q-id + (-L (ni i));g_i )n _i + qsv)o 

= ((L (n , 1) );g_ 1 , n _ 1 + g^)o (2.41) 

Now combining Q, flT3ll) and (ggg) yields 

= -Q-1-1 + ipQ-1,-1 

= -((L( ft) i))+Q_ 1)n _i + Qo^)o + V> / (0Q-l,n-l + ^Qo -1) 

= -((L M )*X*iP) - ((L (nil) )*g_ ltn _i) + V /"(^Q-i.n-i + Wo,-i) (2.42) 
which is one of the desired Hamiltonian flow equations if we identify 

X = — — = d 1 - 1 h d n 



= ^ (2-43) 

o - m 

The identification ( |2.43| ), which is a generalization of the one used for GD hierarchy ( |1.2| ) 
|U, is motivated by the following relation 

= / (t— <fao + ■ • • + - Su n ^ + —5<p + —5*1)) 

J ouq ou n -i ocp dip 

= Tr(Q5C (n>1) ) (2.44) 

where Tr(A) = J tr(A) and tr(A) denotes the ordinary trace of a square matrix. The relation 
( j2.44j ) is actually quite common in the matrix formulation of an integrable hierarchy ||. It 
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is related to the fact that the operation Tr(A) provides a natural scalar product between 
the cotangent space (to which Q belongs) and the tangent space (to which &C( n ,i) belongs) 
of the phase space manifold. In the case of GD hierarchy the identification of this sort can 
be derived from the associated linear system fll8| . However, we have not devised a similar 
proof for the present case. Here, the validity of ( |2.43| ) will be simply justified by the final 
result. 

Now we treat a little differently An equivalent form of fl2.40p is 



which leads to 



= -Q-! + d-\^Q + Q_ 1 , ft _ 1 (L (nil) ) + ) 



Q-i = {Q-i)+ = (9"VQo + d^Q-^L^). 



(2.45) 



(2.46) 



From (|2.37| ) we have 

Q- 1 = [(Qo,-i-X<P)d- 1 }>. n (2.47) 
Now putting ( ggp , (pT47p into we obtain 

£(-k(n,l)) + — {L(n,l)X) + (L( n ,l))+ — {L( n! i)) + (XL( n ,l))+ — [£(n,l)(Qo,-l ~ X(f))d 



-<fi(d 1 ipQ + d 1 Q_i in _iL (ni i ) ) + 
Using (|2.27| ) and (|2.43|) we end up with 



(2.48) 



5H S 



5(L(n,i))+ = (L(„,i)X) + (L( ni i)) + - (L(„,i)) + (II(„,i)) + + [((L( nj i)X) + + L( n>1 )— )9 



-[<f>d- 1 iP(XL {ntl) ) + + 



(n,l)J + 



= [(X(„,i)X) + L (ri) i) - L(„,i)(II(n,i))+] + + [L ( „,i)^<9 V - 09 1 ^L (nil) ]+, (2.49) 

an expected result. 

Finally we turn to ( ^.7| ) and ( |2.13| ) : 



n-l 



iV n _l _i = -Q' n -i _! + 0Q„-i, n _i - - J! u kQk,~l 

k=0 



(2.50) 



Note first that (02|) gives 



Ql,-1 — Qo,-l — 0Qo,n-l 
Q2,-l = Qo-1 — 0Ql,n-l 



0,71-1/ 



Qj,-1 — Qch-1 ~~ 4>Qj-l 



j-l,n-l 



j-2,n-l; 



0,n-l 



(2.51) 



Qn- 



1,-1 



n («-l) 

Vo,-i 



n— 2,n— 1 



n— 3,n— 1 ; 



0,n-l 



,(n-2) 



s 



In particular, 



Q'n-1,-1 — Qo-l ~~ 0Qn-l,n-l — {4>Q n-2,n-l , 



0,n-l, 



v(n-l) 



n— l.ra— 1 



(2.52) 



The operator forms of ( |2.51| ) and ( |2.52|) are, respectively, 

Qj,-i = [&(Q ,-i - X0)] o (j = 0, 1, ■ ■ ■ , n - 1) 

Q'n-1,-1 = [9 n {Q ,-l - X<P)] + <f>Qn-l,n-l 



Putting ( |2.53| ) and ( |2.54|) into ( |2.50| ) we obtain the desired equation : 

n-l 



5(j) 



-[d n (Qo,-i - X<f>)\ Q - J2 Ujl&iQo,-! ~ X(f>)] 

k=0 

[(L M ) + W + f)]„-^f(0^-^) 



Sip- 

{L( n ,l)X<f))o + (L( n> i)— ) - <f) 



Sip 

x , ,5H 5H 



-1,-1 



(2.53) 
(2.54) 



(2.55) 



If we remember identification, 5f = {f,H} mentioned in the paragraph following (|2.14D , 
where / is either of the dynamical variables Ui, <fi and ip, then eqs. ( p.42[ ), ( p.49| ) and (p.55[) 
together define the second Hamiltonian structure of the one-constraint KP hierarchy. In 
conclusion, we have shown that the bracket of Oevel and Strampp || comes out naturally 
from the matrix formulation (which, of course, is basically the AKNS scheme). 



III. GENERALIZATIONS 

Having derived the Hamiltonian structure of the one-constraint KP hierarchy from a 
matrix formulation we now come to the multi-constraint case. The Lax operator for this 
case is defined by 



i=i 



which has matrix representation of the form 



(3.1) 



£(n,m) = Id + U 

( o •■■ o -Vv, 



u 











-1 



\ 



h u Ui 



-1 

• ■ ■ «n-l / 



(3.2) 



(3.3) 
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where U is a (n + m) x (n + m) matrix and the numeration of rows and columns will be 
i, j = —m, — m + 1, • • • , n — 1. Our matrix is connected by the Miura transformation, within 
the generalized Wilson-Drinfeld-Sokolov method, with the matrix used in ref. |TJ]]. We also 

introduce another (n + m) x (n + m) matrix Q such that [Q, £(n,m)] consistent with the form 

of 5C(n t m). We can compute the matrix elements = [Q, £(n,m)]j- an d make it of the form 
<5£(n,m)- The result is the following : 

(1) % < or j < 

AL Q _jg = -Q'_ a _p + 4>/3Q- a , n -l + 1paQo,-/3 (3-4) 

m 

^-a,0 = -Q- a ,0 + «oQ-a,n-l + ^a<5o,0 ~ X! ^PkQ-a-k (3.5) 

fe=l 

iV_ a)i = -<5'_ aji - Q-a,i-l + UiQ-a,n-l + ^aQo,i (z = 1 , 2, • • • , 71 - 1) (3.6) 

iVj _ a = -Qj _ a + (f> a Qi. n -l + Qi+1 ,- a (i =0,1, •••,71 -2) (3.7) 

m n— 1 

^"n-lj-a = ~~ Q'n-l- Q + 4>aQn~l,n-l ~ E ^kQ-k-a — X U kQk-a (3.8) 

fe=l fc=0 

where a, /3 = 1, 2, • • ■ , m. 

(2) i > and j > 

m 

iVi,o = -Qi, + M oQi,n-i + <5i+i,o - X V'fcQi -fc (i = 0, 1, •••,n-2) (3.9) 

fe=i 

m m n— 1 

Nn-lfl = ~Q'n-l,0 + M 0<5ri-l,n-l ~ X/ i'kQn-l-k ~ E 4>kQ-k,0 ~ E U kQk,0 (3.10) 

fc=l fe=l fc=0 

Nij = ~Qi,j - Qi,j-i + UjQ i>n -i + Qi+i,j (*•_?' I' (3-11) 

N n -i f j = —Q'n-l,j ~ Qn-l,j-l + u jQn-l,n-l ~ 
m n—1 

E <^-*y - E Q*J (j = 1, 2, • • • , n - 1) (3.12) 



fc=l k=0 



with 



(3.13) 



N_ jfi = -S^ (3.14) 

N n - 1 - j = 5<j>j j = 1, 2, •••,771 (3.15) 

Nn-ij = <3Mj j = 0, 1, • • • , n - 1 (3.16) 

iVjj = otherwise (3-17) 

The definition of Qi, Q l and iVj are still given by (|2.15| ) and (|2.16|) . Following the steps 
presented in the previous section and imposing the following identifications 

SH 



X = Q 



n-l 



Q-k,n-l — TT (3.18) 
d((>k 

5H 

Qo,-k = —rr (fc = i,2,---,m) 

Oipk 
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we obtain after performing a similar derivation the Hamiltonian structure associated with 
the Lax operator flO| ) : 



&{L(n,m))+ — [{L(n,m,)^)+L( n ,m) ~ -^(n.m) (^-^(n,m))+] + 

+ ^2[L( n ,m)YT~ ^""V* ~ <Pkd~ 1 Y—L( n ,rn)]+ (3.19) 
k=1 OWk 0<Pk 



54>j = {Lu^X^^q + {Li n , m )TT )° (0jTT ^k-r—) (3.20) 

d vi fc=i "^j 

Sipj = -((L M )*X*^j) - {{L (n , m) Y— ) + 51 Tpkj ^ k JT.~^i~^r} ( 3 - 21 ) 



We would like to remark that the Hamiltonian flow equations (|3.19| )- (|3.21|) are, in fact, 
coming from the second GD brackets defined by the Lax operator ( |3.1| ). To see this, let us 
denote tt^- — = si T SH ^ — I - A where A is a differential operator. From the identity 



0L (n;m) ) = tr — — (5 (L(„, m ) + + }^ / ("FT - ^ + "FT* 

£(n,m) d(L (n>m) ) + ^7 5(f) j 6lpj 



we find that the differential operator A satisfies 



(Afc)o = (^)o = |^ (3.23) 



y 3 ^rj 



Now putting L( n m ) and tt = jrj % — h A into the second GD bracket (11.3) and using 

the relation ( |3.23| ), we obtain the Poisson brackets ( p. 19| )-( pT21] ) . 

Note that when we set the next leading coefficient to vanish we have to impose 
the constraint res[-n^ — ,L/ nm )] = 0. Then the resulting bracket becomes the modified GD 

bracket ( |1.5|) defined by L( n ,m) with w n _i = 0. From this modified bracket, it can be shown 
|T5fl that the Hamiltonian structures of the coupled AKNS hierarchy (n = 1, m = 2) and the 
coupled Yajima-Oikawa hierarchy |Tj| (n = 2, m = 2) come out as the special cases. 

In the remaining part of this section, we would like to consider another reduction of the 
KP hierarchy which has Lax operator of the form 

K = 8 n + u^d 71 - 1 + --- + u + a—*— (3.24) 

a — s 



Such system was introduced in the study of matrix model [17] and can be written in the 
form (|2.1[ ) if we equate with each other the integral parts of the two Lax operators : 

00-^ = a—— (3.25) 
a — s 

Then 

a = 0^ s = -(1bV)' (3-26) 

or, equivalently 



11 



ae 



r 



r 



(3.27) 



We can also formulate the Hamiltonian structure associated with K in matrix form by 
introducing 



fd-s -10 
•■■ d -1 








d -1 
■ ■ ■ d + u n -i 



\ 



(3.28) 



and choosing a (n+ 1) x (n + 1) matrix R such that [R, K] consistent with 5K. 
After equating the matrix elements of [R, JC] with the matrix elements of 5JC: 



71- 1) 





= -5s 






StCn-1,-1 


= 5a 








= 5uj 


0' = o,v- 


•,n 




= 


otherwise 





(3.29) 
(3.30) 
(3.31) 
(3.32) 



we can solve Rij in terms of R-x,-i and -Ri,„_i from the constraint equation ( |3.32| ). Thus the 
remaining task is to identify R-±-i and Ri >n -i. However, the following naive identifications 



R-i,-i 



5H_ 
~~5s~ 



R 



i,n—l 



5H_ 

5ui 



0, n-l) 



(3.33) 
(3.34) 
(3.35) 



do not give us the correct Hamiltonian structure. To get correct identifications we note that 
Eq.( p.26| ) can be expressed as a matrix equation: 



where the matrix \1/ is defined by 



(3.36) 



(3.37) 



1/ 



From ( ^.36| ) and the Hamiltonian flow equation 8C(n,i) = [Q, £(n,i)] h is easy to show 
that 5JC = [— \I/ 1 <5\I/ + ^/^Q^/, JC). As a consequence, R is related to Q determined in the 
previous section by 



r = -^~ 1 5^ + ^- 1 g^ 



(3.38) 
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From this relation we find that the correct identification of R-x,-i should be 



R-i, 



OS 



(3.39) 



Substituting ( |3.34| ), ( [3.35D and ( |3.39| ) into the "dynamical" equations, ( |3.29| )-( |3~3"TD , in- 
volving 5a, 5s, and 5L + we obtain the Hamiltonian structure 



5K 4 



(K + X) + K + - K + (XK + ) + + [K + (Xa + + (^)')(d 



-a[{d 



5a 



(3.40) 



5a = [e-f s K + (Xa + + &')ef* s } - a[ef' " S K*JX* + 5 JL) e ~ T % (3.41) 

5a 5s 5a 



(3.42) 



where X = j^-- One can show by explicit calculations that ( |3.40|) -( j3.42|) is indeed the same 



as the second GD bracket defined by K [18 



Rigorously to say, what we obtain here is not really a matrix formulation of the second 
Hamiltonian structure associated with K. This reason is that the identification of _i, 
( |3.39j ), contains the term d~ x 5s which can be determined only after 5s is computed. In 
other words, R is not really the "dual space" of the matrix K, (which contains all dynamical 
variables). The question whether or not a matrix formulation for the second GD bracket 
defined by K remains open. 



IV. EXAMPLES 

In this section, we consider two explicit examples. The first is the cKP hierarchy associ- 
ated with the Lax operator 

£(2,i) =d 2 + u + 0<TV- (4.1) 

It is easy to transform the corresponding eigenvalue equation L^,i)^P = into a matrix form 
by taking 

d -i) \ fcpt 
£ ( 2,i) =[09-1 $= U>2 I (4.2) 



Then £( 2 ,i)$=0 implies 



u d J \(f 3 



iff x - i)ip 2 = (4.3) 
<p' 2 ~ <Ps = (4.4) 
(jxpi + u(f 2 +(p' 3 — (4.5) 



Substitutions of (O) and (O) into (O) gives 
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(<9 2 + u + <j>d~ W2 = (4.6) 

which is precisely L(2,i)¥> = once if = tpi is imposed. 

To consider the Hamiltonian structure associated with £(2,1), let us consider the associ- 
ated Lax equation 

dtC {2 , 1} = [M,£ { 2,i)] (4.7) 

It is necessary to solve the form of M. to make the above equation consistent. 
Writing 

/ a b c \ 
M= d e f (4.8) 
\g h —a — e J 

we have the following explicit expression for the commutator [M, £(2,1)] : 

(ccj) + ipd —aip + cu + ipe —b + tpf 

f<t> + g -dip + fu + h -a-2e 

—(2a + e)0 — ud, —gip — (a + 2e)u — 06, —h — ccj) — fu 

(4.9) 

Setting the matrix (|4.9|) to equal to 

-5ip 

8C {2 ,i) = I I (4.10) 



Su 



we obtain five constraint equations 



—a' + op + ipd = 
-c' - b + ipf = 
-d' + f ( f ) + g = 
-e - dip + fu + h = 

-/'-2e-a = (4.11) 

It is easy to see that a, b, e, g, h can be solved in terms of c, d, f and u, <p, ip 

a = (9 _1 (c0 + ipd) 
b=-c' + -ipf 

9 = d'-f4> (4.12) 
e = -\W + d~\ccP + ipd)) 

h = -\{f" + c<P-d^j)- fu (4.13) 

On the other hand, we have 
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Hence 



5u = —h' — gip — (a + 2e)u — 
50 = —g' — (2a + e)0 — ud 
5ip = b' + (a — e)ip — cu 



Su = (i<9 3 + du + ud)f + i(«90 + 0«9)c - ^(d^j + tpd)d 
50 = ^(90 + 09)/ - ^(0<rV)c - (d 2 + u + 

= -(^ + ^<9)/ -{d 2 + u- |^-V)c + |(^<9" V)d 



If we identify 



/ - — c- — 

5m 50 



(4.14) 
(4.15) 
(4.16) 



(4.17) 
(4.18) 
(4.19) 

(4.20) 



and regard Su, 50 and 8ip as Hamiltonian flows, then we can read off the Poisson brackets 

.1 



{u(x),u(y)} = [^dl + d x u(x) + u(x)d x ]5(x - y) 



{u(x),(f)(y)} 
{u(x),ip(y)} 

{</>(*), <Ky)} 



K x ) 9 x + ~<l>'(x)]6(x - y) 

b{x)d x + -ijj'{x)\8{x - y) 
3 

--(p(x)e(x - y)<p(y) 



{(p(x), ip(y)} = [d 2 x + u{x))S(x -y) + -<p{x)e(x - y)ip(y) 
3 

{ip(x),ip(y)} = --ip(x)e(x - y)ip{y) 



(4.21) 
(4.22) 
(4.23) 
(4.24) 
(4.25) 
(4.26) 



where e(x — y) = d~ 1 5(x — y) is the antisymmetric step function. These are correct brackets 
||, which can be computed from the modified GD bracket (|1.5| ) associated with L = d 2 + 

U + 0<9 _1, 0. 

The second example is to consider the Lax operator of the form 



K = d 2 + u + a 



d — s 



(4.27) 



Even though we would not get a genuine matrix formulation of the second GD bracket 
associated with K as discussed in the previous section. However, we like to check explicitly 
for this simple case that our identifications ( |3.34j) , (|3.35|) and ( |3.39| ) indeed give us a correct 
Hamiltonian structure. 

It is also easy to derive a matrix representation for K 



(4.28) 





fd-s -1 








K = 


d 






02 




\ a u 






Us/ 
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Then from /C$ = we have 



(d - S)0! - 2 = 
02 - 03 = 



ad) i + M0 2 + 







which imply 

[d 2 + udj 2 + a{d- s)- 1 ]dj 2 = 
as desired. The associated Lax equation is given by 

dtK = [M,K] 



(4.29) 
(4.30) 
(4.31) 



(4.32) 



(4.33) 



Writing 



Af 



I ni n 2 n 3 \ 

n 4 tiq 
\ n 7 n 8 n 9 J 



then 



( ari3 + ri4 —n\ + un 3 + sn 2 + 

— STI4 + ariQ + n-j — n 4 + uuq + 11$ 

\ —sn-j + an 9 — an\ — un^ —rij + un 9 — an 2 — un 5 



(4.34) 



-n 2 + sn 3 + n 6 \ 

-n 5 + n 9 
-n 8 — an 3 — uuq J 

(4.35) 



which must be equal to 



-5s 
SIC = ( 

5a 5u 



(4.36) 



There are six constraint equations 



— n 2 — rii + un 3 + sn 2 + n$ = 
-n' 3 -n 2 + sn 3 + n 6 = 
— u'a — sn 4 + an 6 + n 7 = 



—n' 5 — 114 + utiq + 11$ = 



-n 



n 5 + n 9 = 
— n' 9 — 1%$ — an 3 — uuq = 

and three dynamical equations 

5s = n'i — an 3 — n 4 

5a = —n' 7 — sn-j + an 9 — an\ — wn 4 

5u = —n' s — n 7 + un 9 — an 2 — un 5 



(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 



(4.43) 
(4.44) 
(4.45) 
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We can solve n 2 , n±, n 5 , n?, ng, and rig from ( |4.37|) -( p~42|) in terms of rii, n 3 , and uq as follows 



n 2 = —{d - s)n 3 + n 6 

n 4 = (2d(d - s) 2 + 29m - a)n 3 - 2dm - d(3d - 2s)n 6 
™5 = —((9 - s) 2 + u)n 3 + (d - s)n 6 + Tlx 
n 7 = (2(d + s)d(d - s) 2 + 2(d + s)du - (d + s)a))n 3 
-2(3 + s)dn 1 - ((d + s)d(3d - 2s) + a)n 6 
(d(d — s) 2 + du — a)n 3 — (2d 2 — ds + u)n 6 — drii 



n 9 



-((d - s) 2 + u)n 3 + (2d - s)n 6 + rii 



Now substituting ( |4.46| )-( ^oTD into the dynamical equations ( |4.43| )-( |03| ) and 
identifications 



6H ^x-ix 

Tlx = h O x OS, 

OS 



n 3 



5H_ 

5a 



5H_ 

5u 



(4.46) 
(4.47) 
(4.48) 

(4.49) 
(4.50) 
(4.51) 

using the 
(4.52) 



we obtain the expected Poisson brackets 18 



{«( 


x) 


,u(y)} = 


+ 2u(x)d x + u'(x)]5(x - y) 


(4.53) 


{u 


x) 


,a(y)} = 


[3a(x)d x + 2a'(x)]5(x — y) 


(4.54) 


{u 


(x) 


Xv)} = 


\h 2 x + s (x)d x )8(x - y) 


(4.55) 


{a 


x) 


,a(y)} = 


\(Aa(x)s(x) + 2a'(x))d x + a"(x) + 2(a(x)s(x))']5(x - y) 


(4.56) 


{a 


(x) 


Xv)} = 


\(d x + s(x)) 2 d x + u(x)d x ]5(x - y) 


(4.57) 


{s 


(x) 


Xv)} = 


-d x 5(x-y) 


(4.58) 



V. CONCLUDING REMARKS 

We have studied the Hamiltonian structure of the cKP hierarchy in the matrix formu- 
lation. We showed that the procedure for the GD hierarchy can be gone through without 
difficulty for the cKP hierarchy. We have not only reproduced Oevel and Strampp's result 
for the one-constraint case but also generalized it to the multi-constraint case. We further 
showed that this Hamiltonian structure is nothing but the GD bracket. Hence our result 
again confirms the fact that the GD bracket can be properly restricted to the Lax operator of 
the form ( |3.1| ) . We have also discussed the matrix formulation of the Hamiltonian structure 
of the one-constraint KP hierarchy in the form given by (|3.24j ). However, we did not obtain 
a genuine matrix formulation in the usual sense. Nevertheless, we found that one can still 
compute the Hamiltonian structure from our matrix formulation once a proper modification 
is made. 

Finally, we want to remark that the constrained modified KP hierarchy can be ob- 
tained from the Lax operator Q2.ip via the gauge transformation 

K = _1 L( n) i)0 

= d n + +---+Vo + d- l v^ (5.1) 
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which satisfies the hierarchy equation 



f = (5-2) 



The bi-Hamiltonian structure associated with ( |5.2| ) has been obtained M from the lifted 
bracket of the cKP hierarchy by gauge transformation. It would be interesting to know 
whether or not a matrix formulation of this bi-Hamiltonian structure exists. Following the 
spirit of the present work we have worked out this structure in matrix formulation for n = 2 
and 3, but a general construction is yet to be given. We hope to report the results in this 
direction in the near future. 
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